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We present a systematic study of the phase-coherent dynamics of a superconducting three- 
Josephson-j unction flux qubit. The qubit state is detected with the integrated-pulse method, which 
is a variant of the pulsed switching DC SQUID method. In this scheme the DC SQUID bias current 
pulse is applied via a capacitor instead of a resistor, giving rise to a narrow band-pass instead of 
a pure low-pass filter configuration of the electromagnetic environment. Measuring one and the 
same qubit with both setups allows a direct comparison. With the capacitive method about four 
times faster switching pulses and an increased visibility are achieved. Furthermore, the deliberate 
engineering of the electromagnetic environment, which minimizes the noise due to the bias circuit, 
is facilitated. Right at the degeneracy point the qubit coherence is limited by energy relaxation. We 
find two main noise contributions. White noise is limiting the energy relaxation and contributing 
to the dephasing far from the degeneracy point. 1//- noise is the dominant source of dephasing in 
the direct vicinity of the optimal point. The influence of l//-noise is also supported by non-random 
beatings in the Ramsey and spin echo decay traces. Numeric simulations of a coupled qubit-oscillator 
system indicate that these beatings are due to the resonant interaction of the qubit with at least 
one point-like fluctuator, coupled especially strongly to the qubit. 



; I. INTRODUCTION 

In the past two decades it has become evident that 

\ quantum mechanics can give rise to new fascinating pos- 
sibilities for communication and information processing. 
Based on quantum two-level systems (qubits) instead of 
conventional classical bits, secure quantum communica- 
tion protocols and quantum algorithms have been pro- 
posed. The latter promise a significant speed-up of cer- 
tain computational tasks, such as the factorization of 

\ large numbers, by making use of a massive quantum 
parallelism and entanglementii^i^i^i^i. Subsequently, 
there have been intensive research efforts aiming at the 
development of adequate hardware concepts^. The suc- 
cessful implementation of qubits has first been proposed 
and demonstrated for microscopic systems (e.g., NMRi^, 
trapped ionsiSiii, and cavity-QED systemsi^iHiiiii^ii^) . 

\ Although these systems possess sufficiently long decoher- 
ence times, they have drawbacks regarding their scalabil- 

' ity to large architectures, which are required for practi- 
cal quantum computing. In contrast to that, scalability 
is a specific advantage of solid-state systems. Therefore, 
solid-state quantum circuits have attracted increasing in- 
terest giving rise to a rich field of theoretical and exper- 
imental investigations in the context of quantum infor- 
mation processing. The fabrication process of solid-state 
quantum circuits is strongly facilitated by the use of well 
established techniques from micro- and nanoelectronics 
such as lithography and thin-film technology. Obviously, 
this results in a high degree of tunability and potential 
scalability to larger units. Furthermore, quantum sys- 



tems clearly have a great potential in solid-state electron- 
ics because the ongoing miniaturization of integrated cir- 
cuits requires the deliberate use of quantum-mechanical 
effects in the near future. 

Among solid-state systems, superconducting devices 
are particularly interesting. Quantum information pro- 
cessing based on superconducting circuits exploits the in- 
trinsic coherence of the superconducting state. In fact, 
this state is separated by an energy gap from the nor- 
mal conducting states and, therefore, relatively long co- 
herence times can be achieved. In general, supercon- 
ducting quantum circuits consist of a single or multiply- 
connected superconducting lines intersected by Joseph- 
son junctions. Quantum information can be stored in 
the number of superconducting Cooper pairs (e.g., the 
charge qubi t^^i^^'^^ and the quantroniu m^°i^^ ) , in the di- 
rection of a circulating persistent current (e.g., the three- 
junction flux qubi1,22i23^24^^^ ^^^^ ^Yve RF-SQUID- 

based flux qubili^) or in oscillatory states (e.g., the phase 
qubil!2^i^i^>^i^i2ii2^). To a large extent, superconduct- 
ing qubits are protected from the unwanted interaction 
with environmental degrees of freedom by the super- 
conducting gap. Furthermore, design-dependent internal 
symmetries can reduce the influence of noise arising from 
the control and readout circuitry. Nevertheless, the un- 
controlled loss of coherence still represents a key issue. In 
the most simple theoretical description^i^i decoherence 
in a quantum two-level system is described in terms of 
two rates or times. The first is the longitudinal or energy 
relaxation rate Fi = Tj^^, which describes the transi- 
tions between the two qubit states due to high-frequency 
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fluctuations. The second is the transverse or dephasing 
rate = T^^, which describes the loss of phase coher- 
ence within the same qubit state caused by low-frequency 
noise. While Fi can be measured directly, has to 
be extracted from the experimentally accessible quantity 
F2 =Fi/2 + F^. 

The coherence time of a quantum state is one of 
the key figures of merit for practical quantum infor- 
mation devices. In fact, it determines the number of 
qubit operations that can be performed without errors. 
The presently achieved coherence times in superconduct- 
ing qubits (ranging from approximately 10 ns to 5 fis) 
are not sufficient for the realization of more complex 
quantum circuits. For this reason, decoherence due 
to the coupling of the qubit to the environmental de- 
grees of freedom or due to noise introduced by the con- 
trol and readout circuitry has been addressed in several 
recent experimenta P^i^^'^''i'^^'^^i'^^'^°i^^''^^ and theoretical 
studiea^ii^i^i^i^. However, a thorough understanding 
of the decoherence processes and detailed knowledge of 
the origin of the noise sources is still missing, making 
the further clarification of possible sources of decoher- 
ence highly desirable. In this sense, a promising path 
to follow is to study the dependence of the character- 
istic decay times (energy relaxation and dephasing) on 
external parameters such as the applied magnetic flux in 
superconducting flux qubit a^^i^^'^^ . 

It is noteworthy to mention that the interaction of su- 
perconducting qubits with microwave resonators has re- 
cently attracted increasing attention. It has turned out 
that the qubit-resonator interaction is the circuit equiv- 
alent of the atom-photon interaction in optical cavity 
QEDi^iHiiiii^i^. In other words, the formalism of cavity 
QED has been successfully transferred to the realm of su- 
perconducting systems, which can now be treated with 
the standard tools developed in quantum optics. This 
new field is referred to as circuit QED and has given 
rise to a set of important method for the manipulation 
and dispersive readout of qubit cj'^''''^^i^^'^°i^"'^'^^ . Partic- 
ularly interesting for the purposes of the present work 
are the investigations on coherent dynamics presented in 
Refs. [5^ and [SJ. Together with our work they allow for 
further insight in the issue of nearly-resonant spurious 
resonators or two-level systems interacting with super- 
conducting qubits^. 

In this paper, we present a detailed study of both en- 
ergy relaxation and dephasing rates of a superconduct- 
ing three- Josephson-junction flux qubit as a function of 
the applied magnetic flux close to the degeneracy point. 
For this type of qubit several different readout methods 
have been proposed and successfully implemented. They 
range from the simple switching-DC-SQUID method^ 
to more sophisticated techniques such as the inductive 
readoui^^i^i^i^ or the bifurcation amphfie r^^i^° i^i^. 
The two latter have the potential to achieve quantum 
non-demolition measurements^ and very large signal vis- 
ibility, which are very important features for future ap- 
plications. In this study we have chosen a switching 



DC SQUID readout, which enables us to study the ef- 
fect of two fundamentally different electromagnetic envi- 
ronments (capacitive and resistive) on the decoherence 
of one and the same flux qubit easily. Furthermore, 
the switching-DC-SQUID readout is attractive because 
of its simple technical implementation (e.g., no cold am- 
pliflcrs or other sophisticated high-frequency circuits are 
needed). In our measurements of the coherent dynamics 
of a superconducting flux qubit we apply the bias pulse to 
the DC SQUID detector via a capacitor (C-Bias). Such 
a setup reduces low-frequency bias current noise. We 
experimentally demonstrate that the decoherence of our 
qubit at the optimal point is dominated by relaxation, 
i.e., by high-frequency noise. Nevertheless, we still find 
the signature of l//-noise in the vicinity of the optimal 
point. We compare our results to those obtained from 
resistive-bias (i?-Bias) measurements on the same qubit 
and find a similar 1// spectral noise density, despite the 
different bias line filtering conditions. This demonstrates 
that low-frequency noise from the bias line is not the 
dominating noise source in our system. The relaxation 
time at the degeneracy point is slightly reduced for the 
capacitive method. We are able to exclude the different 
bias line configuration as the reason for this Ti-time re- 
duction. Instead, we attribute it to a change in the high- 
frequency setup when switching between the two meth- 
ods. Furthermore, we show that the resonant interaction 
of the investigated flux qubit with two-level fluctuators 
is responsible for non-random beatings in the measured 
Ramsey and spin echo traces. Together with the observed 
weak influence of low-frequency bias noise on the qubit 
coherence, this suggests that an ensemble of fluctuators 
on the sample chip and close to the qubit is the main 
source of the measured l//-noise. Additionally, we find 
that in our setup, compared to the resistive-bias measure- 
ments, the capacitive-bias method allows for four times 
faster switching pulses and helps to obtain an increased 
visibility. 

The paper is organized as follows. In Sec. |TT1 we first 
introduce the fabrication process and the basic properties 
of the investigated three-junction flux qubit, the exper- 
imental setup, and the general measurement concepts. 
Then, in Sec. lIIIl we describe the details of the capacitive- 
bias readout scheme, which is based on the switching 
of a DC SQUID. Here, in particular the differences to 
the usual resistive-bias setup are discussed. In Sec. IIV[ 
we show the results of the capacitive-bias measurements 
on the phase coherent dynamics of our flux qubit. We 
present a detailed quantitative analysis of the results of 
these measurements, which have been performed both 
in the frequency (microwave spectroscopy) and in the 
time domain (energy relaxation, Ramsey fringes, and 
spin echo decay). We discuss the effect of the different 
electromagnetic environments by way of comparing the 
capacitive-bias results to those obtained from resistive- 
bias measurements on the same qubit. In this way, we 
avoid possible ambiguities arising from a comparison of 
different samples due to the spread in fabrication and/or 
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the rich variety of different decoherence sources. From 
the measured flux dependence of the energy relaxation 
and dcphasing rates we are able to quantify the contri- 
butions of 1//- and white noise to the total flux noise 
spectral density affecting the qubit. Finally, in Scc.|Vl we 
discuss the origin of non-random beatings in the Ramsey 
and spin echo decay traces. We show that these beatings 
cannot be caused by imperfect control pulses or undesired 
probe frequency detuning because in our experiments we 
apply the so-called phase-cycling method^l. Comparison 
with numerical simulations strongly suggests that these 
beatings originate from the interaction of the flux qubit 
with one (or a few) point-like harmonic oscillators or two- 
level fluctuators. 



II. THE SUPERCONDUCTING QUANTUM 
CIRCUIT 

In this section, we introduce the working principle and 
characteristic parameters of our work-horse, the super- 
conducting three-Josephson-junction flux qubit. We in- 
troduce the pulsed switching DC SQUID readout tech- 
nique, placing special emphasis on the capacitive-bias 
method. Furthermore, the control of the qubit state with 
microwave pulses as well as the filtering of the measure- 
ment lines are discussed. 



A. The superconducting three-Josephson-junction 
flux qubit 

The investigated sample is fabricated on a 3.5 x 
3.5 mm^ Si02-passivatcd Si substrate. Standard electron 
beam lithography and aluminum thin-film technology are 
used. The sample chip is surrounded by a T-shaped 
printed circuit board (PCB). Both sample chip and PCB 
are placed inside a gold-plated copper box, which is con- 
nected to the mixing chamber stage of a dilution refrig- 
erator. All experiments are performed at the base tem- 
perature of the fridge, which is approximately 50 mK. In 
the following, we briefly describe the main components of 
the sample and summarize their characteristic parame- 
ters and function. A sketch of the sample layout is shown 
in Fig. m 

The three-Josephson-junction flux qubit^^ consists of 
a square-shaped superconducting aluminum loop inter- 
rupted by three nano scale Al/AlO^j/Al Joscphson junc- 
tions. The area of two of these junctions is chosen to 
be the same (0.03 /im^), whereas the third one is de- 
signed to be smaller by a factor adosign = 0.7. The 
qubit (and DC SQUID) junctions are fabricated using the 
shadow evaporation technique. The Josephson junctions 
are characterized by their Josephson coupling energy 
Ej = /c$o/27r and their charging energy Eq = e^/2Cj. 
Here, $o is the flux quantum, e the electron charge, and 
Ic and Cj are the critical current and capacitance of the 
junction, respectively. The specific capacitance of our 



junctions was determined to be^ Cs = 100 ± 25 fF / fim^ . 
When the ratio Ej/Ec is approximately 50, the device 
is expected to behave as an effective quantum two-level 
system (qubit). 

The effective Hamiltonian of a flux qubit can be writ- 
ten as22d3i65 



^qb — — o — ^ 



A. 



(1) 



where and ax are the Pauli matrices. For van- 
ishing tunnel coupling (A = 0) the two qubit states 
correspond to the classical states |R) and |L) with 
clockwise and counterclockwise persistent currents Ip = 
±Ic\/T^-j2a)~^ circulating in the loop2^. They are sep- 
arated by the flux-dependent energy e($x) = 2/p(5<I'x, 
where S^x = — ^(n) is the applied magnetic flux 



bias relative to the degeneracy point $(„) = 



and n is an integer, 
erated close to $x 



Our sample is designed to be op- 
$(1). For finite coupling (A > 0), 
we obtain superpositions of |R) and |L). This results 
in new qubit eigcnstates |0) and |1) whose energy dif- 
ference Eoi = ^e($x)2 + A2 = ^(2/p(5$x)2 -I- has 
a hyperbolic flux dependence. At the degeneracy point 
((5$x = 0, e($x) = 0) the qubit is protected from de- 
phasing because Eqi is stationary with respect to small 
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FIG. 1: Sketch of the sample layout. The Al/AlOx/Al 
Josephson junctions of qubit and DC SQUID are represented 
by the symbol Kl. Vin and Vout are the bias pulse and the 
switching signal of the DC SQUID, respectively . Csh is the 
effective shunting capacitance of the DC SQUID. The boxes 
titled "bias" and "read" represent the circuit elements used to 
engineer the electromagnetic environment of the qubit. They 
are located partially on-chip and partially on the PCB sur- 
rounding the sample chip, "microwave & ASP" denote the 
microwave control pulse pattern and the adiabatic shift pulse, 
which are attenuated at low temperature and coupled to the 
qubit via an on-chip antenna (coiled shape). 
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variations of the control parameter 5$x- Therefore, this 
point represents the optimal point for the coherent ma- 
nipulation of the qubit. Also, the qubit eigenstate at 
the degeneracy point is an equal superposition of |L) 
and |R), i.e., the expectation value of the persistent 
current vanishes. Far away from the degeneracy point 
(e(*i'x) 3> A) the qubit behaves as a classical two-level 
system. This example clearly shows the flexibility of- 
fered by solid-state-based qubits due to their high degree 
of tunability. As detailed in Sec. HVl for our device we 
find A//i ~ 4 GHz, Ej/Eq ~ 50, and a critical current 
density Jc ~ 1300 A/cm^. This means that at an oper- 
ation temperature T ~ 50 niK the condition fceT <^ A 
required for the observation of quantum effects is well 
satisfied. 

When describing the influence of fluctuations 6lj on the 
qubit (cf. Sec. lIVp . it is convenient to express the Hamil- 
tonian of Eq. ([T|) in a two-dimensional Bloch vector repre- 
sentation: i/qb = hu:&/2. Here, <t = {a±,a\\) = {^x,^z) 
and u> = (ti>^,W||) = (a, e($x)) ■ The representa- 
tion of the Bloch vector lj in the qubit energy eigen- 
basis is obtained by multiplying u) with the rotation 

matrix D = | ^'^^/i^ I from the left. The Bloch 

Vsmy cost^ J 

I 



angle 9 is defined^ via tan9 = A/e($x)- This re- 
sults in sm9 = A/vqi and cos 9 ~ e(<i>x)/foi, where 
i/Qi = cjoi/Stt = Eoi/h = y/A'^ + e(<i>x)2//i is the qubit 
transition frequency. At the degeneracy point the angu- 
lar transition frequency becomes wa = A/S. 



B. Manipulation of the qubit state 

The qubit control is achieved by varying the con- 
trol parameter (5$x and, simultaneously, applying a suit- 
able sequence of microwave pulses Acoa{2TTi't + (j)) with 
frequency ~ i^oi- A single microwave pulse results 
in a rotation of the qub it state vector by an angle 
51 — 2-KVYit\/\ -I- (i5/;/r)^ on the Bloch sphere, where 
t is the pulse duration and the Rabi frequency v^i = 
vyiXA) is a function of the pulse amplitude A. The 
relative phase <\) of the pulse and the detuning b = 
V — i/Qi determine the rotation axis v = (wi,U2,i'3) = 
(z/R cos <j>, sin (j>, S) I + S^. Mathematically, we can 
describe this rotation with the matrix 



cos n + vf(i — COS ri) 

V2V1 (1 — cos O) + V3 sin f2 
(1 — cos O) — V2 sin f2 



f 1112(1 — cos 51) — 

cos 51 + V2{1 - 
V3V2{1 — cos 51) 



U3 sin 51 
cos 51) 
- vi sin 51 



ViV3{l — cos 51) + V2 sin 51 
V2V!i{l — cos 51) — vi sin 51 
cos 51 -f w|( 1 — cos 51) 



(2) 



For the case 5 ~ {) and = 0, Eq. ([2]) describes a rotation 
by an angle i-Kv^t about the x-axis. When introducing a 
finite relative phase cj) the orientation of the rotation axis 
within the x, y-plane changes. Finite detuning results in a 
change of the rotation angle and a tilt of the rotation axis 
out of the X, y-plane. In the absence of any microwave 
radiation [vii = 0) the qubit evolves freely, i.e., its state 
vector rotates about the z-axis of the Bloch sphere with 
a frequency 51froc = ZttJ • t (since we work in a frame 
rotating with frequency v). The corresponding rotation 
matrix _Rr(flfroo) is obtained from Eq. ([2]) by choosing 
v = v, ="(0,0,1). 

The microwave pulses are applied to the qubit via an 
on-chip antenna, which is implemented as a coplanar 
waveguide transmission line shorted at one end. The 
qubit is coupled inductively to this antenna. From a 
FastHenry^ simulation the mutual inductance between 
the antenna and the qubit is determined to be iWmw.qb — 
73 fH. The applied microwave radiation is cooled by 
means of a 10 dB and a 3 dB attenuator, which are ther- 
mally anchored at a temperature of 4K and 100 mK, 
respectively (cf. Fig. [1]). The initial qubit state, which 
in our measurements is always the ground state, is pre- 
pared by waiting for approximately 300 /iS. This time 



is much longer than the energy relaxation time Ti of 
the system. For the frequency and time domain experi- 
ments, the pulse sequences shown graphically in Fig.[2jb) 
are used for the qubit manipulation. For the measure- 
ments in the frequency domain (microwave spectroscopy) 
the qubit is saturated to an equilibrium mixed state by 
means of a sufficiently long microwave pulse followed by 
the readout. In the time domain, driven Rabi oscilla- 
tions are recorded by measuring the qubit state as a 
function of the duration of a single microwave pulse of 
fixed amplitude A. From the measured Rabi frequency 
I'R of these oscillations we determine the durations 
and of tliG TT- and 7r/2-pulses, which rotate the qubit 
state vector by the corresponding angles. In our ex- 
periments we typically choose the microwave pulse am- 
plitude such that = ~ l/Si'R — 8 ns. Then, 
the energy relaxation time is determined by exciting the 
qubit with a 7r-pulse and, subsequently, recording the 
decay of {dz) as a function of the waiting time. Fi- 
nally, the evolution of {^x,y) is probed by the sequence 
7r/2-pulse-wait-7r/2-pulse-readout (Ramsey experiment) 
or the sequence 7r/2-pulse-wait-7r-pulse-wait-7r/2-pulse- 
readout (spin echo experiment). For the latter sequence 
low-frequency phase fluctuations are canceled because of 
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the rcfocusing effect of the intermediate 7r-pulse. As dis- 
cussed in more detail below, certain beatings of the spin 
echo signal can be corrected by recording several traces 
in which the control pulses have different relative phases 



C. Readout of the qubit state 

In order to read out its state, the qubit is surrounded 
by a slightly larger square-shaped aluminum loop con- 
taining two Al/AlOjj/Al Josephson junctions, a so-called 
DC SQUID (cf. Fig. [J). The DC SQUID is sensitive to 
the Hux difference generated by the persistent currents 
flowing in the qubit loop. In our design, the DC SQUID 
is coupled to the qubit via a purely geometric mutual 
inductance Msq^qb = 6.7 pH. Differently from other flux 
qubit dcsign n^^i^^ , there is no galvanic connection be- 
tween DC SQUID and qubit in our sample. This is ex- 
pected to reduce the effect of asymmetry-related issues 
as well as the detector back-action on the qubit. In fact, 
we do not find any measurable bias current dependence 
of the qubit decay time, as recently recently reported for 
shared-edge dcsign a^^i*^^ . The lines used for biasing and 
reading out the DC SQUID detector are heavily filtered 
against noise in the gigahertz range using a combina- 
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FIG. 2: (Color online) (a) Adiabatic shift pulse (ASP) read- 
out, (b) Microwave control pulse patterns for the experiments 
discussed in this work. The boxed values denote either the 
pulse duration t or the corresponding rotation angle of the 
qubit state vector on the Bloch sphere. Free evolution times 
are denoted by the symbol r. In the multi-pulse sequences, 
(j>i and (f>2 are the pulse phases relative to the initial pulse 
necessary for the phase-cycling technique. 



tion of copper powder filters and stainless steel ultra-thin 
coaxial cables at the mixing chamber temperature level. 
The details of the filtering of the measurement lines are 
different for resistive- and capacitive-bias schemes and 
are explained in detail in Sec. IIIII 

The detection of the qubit state (|0) or |1)) is straight- 
forward far away from the degeneracy point. In this re- 
gion the energy eigenstates coincide with the persistent 
current states |L) and |R). The flux generated by these 
currents is then detected with the DC SQUID. However, 
in the vicinity of the degeneracy point the qubit eigen- 
states |0) and |1) cannot be distinguished in this way 
because they become nearly equal superpositions of \L) 
and \R). In other words, the expectation value {IpCfz) 
of the persistent current circulating in the qubit loop 
and thus the flux signal vanishes. To circumvent this 
problem we employ the adiabatic shift pulse method [cf. 
Fig. [2la)] . This method exploits the possibility to per- 
form the readout process sufficiently far away from the 
degeneracy point by applying a control pulse, which adi- 
abatically shifts the qubit in and out of the region around 
the degeneracy point. In contrast to the quasi-static flux 
bias at the readout point, which is generated by a super- 
conducting coil located in the helium bath of our dilu- 
tion refrigerator, the shift pulse is applied to the qubit 
through the on-chip microwave antenna (cf. Fig.[T]). The 
total control sequence (cf. Fig. ^ for initialization, ma- 
nipulation, and readout of the qubit can be summarized 
as follows: First, the qubit is initialized in the ground 
state at the readout point far away from the degener- 
acy point by waiting for a sufficiently long time. Then, 
a rectangular (0.8 ns rise time) adiabatic shift pulse to- 
gether with the microwave control sequence is applied to 
the qubit via the on-chip antenna. In this way the qubit 
is adiabatically shifted to its operation point, where the 
desired operation is performed by means of a suitably 
chosen microwave pulse sequence. Finally, immediately 
after end of the microwave pulse sequence, the qubit is 
adiabatically shifted back to the readout point, where 
the state detection is performed. Note that, in order to 
avoid qubit state transitions, the rise and fall times of 
the shift pulse have to be long enough to fulfill the adia- 
batic condition but also short enough to avoid unwanted 
relaxation processes. 

In the readout process performed right after shifting 
the qubit back to the readout point, the DC SQUID 
is biased with a current pulse of an amplitude just be- 
tween the two switching currents corresponding to the 
qubit states. Depending on the actual qubit state, the 
DC SQUID detector either remains in the zero-voltage 
state or switches to the running-phase state. Only in the 
latter case a voltage response pulse can be detected from 
the readout line. The DC SQUID voltage signal is ampli- 
fied by means of a room temperature differential amplifier 
with an input impedance of 1 M£l against a cold ground 
taken from the mixing chamber temperature level. 

The DC SQUID is shunted with an Al/AlO^/Al on- 
chip capacitance Csh = 6.3 ± 0.5 pF. This capacitance 
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also behaves as a filter and, in combination with the other 
biasing elements (either a resistor or a resistor-capacitor 
combination), creates the electromagnetic environment 
of the qubit. Previous studies^ have shown that a purely 
capacitive shunt results in a much smaller low-frequency 
noise spectral density compared to an i?C-type of shunt. 
This is crucial since one has to reduce the environmen- 
tal noise as much as possible in order not to deteriorate 
the qubit coherence times. The resistive part of the bias- 
ing circuit helps to damp modes formed by the shunted 
DC SQUID and the parasitic inductance/capacitance of 
its leads. Ideally, it should be placed as close as possible 
to the shunted DC SQUID. 

Finally, after a typical "single-shot" measurement se- 
quence as described above, the response signal of the 
DC SQUID is binary (zero or finite voltage state), de- 
pending on the qubit state. In the experiments we mea- 
sure the switching probability Pgwj which is the aver- 
age over several thousands of single-shot measurements. 
For a proper bias current pulse height the qubit state is 
encoded in the value of the switching probability. In 
the best case, the ground state would correspond to 
= 0% and the excited state to P™'"'' = 100%, or 
vice versa. In reality, however, due to noise issues the ac- 
tual visibility P^^^^ — -Pg™'" is usually significantly smaller 
than 100%. 



III. CAPACITIVE BIAS: THE INTEGRATED 
PULSE METHOD 

In this section, we present our qubit readout scheme, 
which is based on a capacitive instead of the standard 
resistive bias for the DC SQUID. We also refer to this 
method as the integrated pulse method because the volt- 
age pulse sent to the DC SQUID detector is the time 
integral of the desired current bias pulse. Furthermore, 
the consequences arising for the filtering of the DC lines 
are discussed. We note that a similar readout method, 
also based on a capacitive bias, has been proposed in 
an effort to design and implement a switching detector 
for Cooper-pair transistor and Quantronium circuit a^^'^^ . 
However, there are remarkable differences to the work 
presented here. Firstly, because of the on-chip shunt- 
ing capacitor we do not encounter the complication of 
frequency-dependent damping. Even for low impedances 
the readout DC SQUID in our experiments rcmaines un- 
derdamped (DC SQUID quality factor Q ~ 4 at approx- 
imately 50 f2). Secondly, we are able to measure the co- 
herent dynamics of our flux qubit (cf. Sec. IIV[) . 

The readout of the qubit is performed with a 
switch&hold measurement technique, which makes use 
of the hysteretic current-voltage characteristic of the 
DC SQUID detector. In the resistive-bias detection 
scheme, a voltage pulse is generated at room tempera- 
ture, fed into the DC SQUID bias line, and then trans- 
formed into a current pulse via a cold 1.25 kJl bias resis- 
tor [cf. Fig. [31[a)]. According to Ohm's law, the voltage 



pulse must have the same shape as the desired bias cur- 
rent pulse. The latter is composed of a short rectangular 
switching pulse (of duration ~ 60 ns) immediately fol- 
lowed by a much longer hold pulse (~ 1 fj,s) of smaller 
amplitude, as shown in Fig. [3l^a). As already explained 
in Sec. [TTl the switching pulse height is chosen such that 
the DC SQUID either switches to the voltage state or 
stays in the zero voltage state, depending on the state 
of the qubit. Hence, the length of the switching pulse 
determines the time resolution of the switching event de- 
tection. The hold pulse level is chosen to be just above 
the value of the DC SQUID retrapping current. Conse- 
quently, if the DC SQUID switches into the voltage state 
during the switching pulse it will not switch back to the 
zero voltage state during the hold pulse. In this way, 
the voltage is sustained for a sufficiently long time in- 
terval allowing the use of a room temperature amplifier 
with reduced bandwidth. At the same time quasiparticle 
generation is minimized. 

In the resistive-bias experiments the bias lines of the 
readout DC SQUID arc severely low-pass filtered [cf. 
Fig. [3ja)] . This is required to reduce the high-frequency 
noise spectral density Sui{uJoi = Eqi/H), which is respon- 
sible for the energy relaxation. In addition, filtering is 
necessary to eliminate part of the low-frequency noise. 
In this manuscript, the term low-frequency noise refers 
to noise which has a large spectral density at frequencies 
much smaller than the qubit level splitting wa — 4 GHz. 
The dephasing time of the qubit is mainly determined 
by this low-frequency environmental noise spectral den- 
sity Sui {uj ^ 0). Consequently, the cutoff frequency of the 
low-pass filter should be chosen as low as possible in or- 
der to attenuate low-frequency noise strongly. However, 
a lower limit is set by the requirement that the readout 
pulse has to be sufficiently short to avoid any deteriora- 
tion of visibility. For this reason, the filter cutoff cannot 
be chosen smaller than about 10 MHz in the resistive- 
bias measurements. In other words, in the resistive-bias 
setup the very low-frequency noise below 10 MHz passes 
unaffected to the sample. 

The integrated pulse readout reduces the effect of low- 
frequency noise from the DC SQUID control lines, es- 
pecially from the bias current line. In a qubit limited 
by this type of noise, the use of a capacitive bias should 
improve the dephasing time as compared to the resistive- 
bias case. The reason is that a band-pass filter instead 
of a low-pass filter configuration is achieved by replac- 
ing the DC SQUID bias resistor with a 0.5 pF capacitor. 
Consequently, only a narrow frequency band (enough to 
permit sufficiently fast readout pulses) is allowed to pass, 
whereas noise at lower and higher frequencies is strongly 
suppressed. In this way, the rise time of the switching 
pulses applied to the readout DC SQUID is reduced by a 
factor of six (i?-bias: 60 ns, C-bias: 10 ns). In turn, this 
should result in a better time resolution and an increased 
visibility. In actual experiments, we observe an improve- 
ment from 20-25% for the resistive-bias setup to 28-36% 
for the capacitive-bias configuration (cf. inset of Fig.[n|). 
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FIG. 3: (a) Top: Resistive-bias setup. "LPFIO. 7(250)" denotes a commercial low-pass filter with 10.7(250) MHz cutoff frequency. 
"SS/CPF" represents Im of ultra-thin (00.33mm) stainless-steel coaxial cable followed by a copper powder filter (bias line: 
25 cm wire length; readout line: 100 cm wire length). Qubit, DC SQUID, and shunting capacitor Csh are indicated with 
a single cross (x). Solid and broken lines represent high-bandwidth semirigid 01.2 mm CuNi/Nb and narrower bandwidth 
stainless-steel braided flexible coaxial cables, respectively. The bias voltage pulse is attenuated by 40 dB at 4K. The bias 
resistors are -Rabias = 250^-1- lkf2 and -Rarcad = 2.25kn-f 3kf2 (on-chip -I- off-chip) . Bottom: Switch&hold readout pulses for 
the resistive-bias setup. Note that 60ns are the width of the portion of switching pulse which exceeds the hold level. Only there 
switching events can be induced, (b) Top: Capacitive-bias setup: The cutoff frequency of the room temperature bias line filter 
is reduced to 100 MHz. The voltage pulses require a weaker attenuation of 3 dB at 4 K. The series capacitors Cbias = 0.5 pF and 
Cread = 470 pF are inserted into the bias and readout lines, forming a band-pass filter. Small series resistors (i?cbias = 511 f2 
and 7?croad = 1.5 kS7) still provide sufficient damping of parasitic external modes. Aside from Csh, all resistors and capacitors 
are placed off-chip, on the PCB. Bottom: Switch&hold readout pulses for the integrated-pulse setup: The voltage pulse (dashed 
line) is the time integral of the desired current pulse (solid line). An approximately 300 /is long pullback section is required. In 
the actual voltage pulses (dotted line) another kink is introduced to avoid discharging effects of the capacitors. 



Since the incoming voltage pulse is differentiated by the 
bias capacitor, the time integral of the desired current 
pulse has to be applied as shown in Fig. E^b). 

In Fig. m we compare the pulse shapes and filter char- 
acteristics of the two different bias schemes. Evidently, 
the integrated pulse readout provides a band-pass filter- 
ing of the bias line and allows for approximately four 
times shorter switching pulses than the resistive-bias 
setup. In part, this is achieved by replacing the 10.7 MHz 
commercial low-pass filter utilized for the resistive-bias 
readout by a 100 MHz fiher [cf. Fig.E^b)]. Note that for 
the integrated pulse setup the different levels of switching 
and hold part of the current pulse correspond to differ- 
ent slopes of the voltage pulse in these sections. Since 
the readout sequence has to be repeated many times, 
the voltage pulse amplitude has to be reduced back to 
zero after the end of the hold section. This so-called 
pullback section has to be much longer (~ 300 /is in our 
experiments) than the hold section in order to avoid any 
switching of the DC SQUID in the negative current direc- 
tion. We note that in the experimental implementation 
the amplitude of the differentiated current pulse falls off 
the switching to the hold level within a time constant of 
a few tens of nanoseconds. This is caused by the dis- 
charging of the shunting capacitor Csh and can be com- 
pensated by introducing a kink in the integrated voltage 
pulse [cf. Fig. [3^b)]. Moreover, we also place a capacitor 



Croad in the voltage output line. Since small capacitors 
strongly attenuate the outgoing signal, the value of Cread 
is limited to 470 pF. 

The microwave control signals applied to the sample 
via the on-chip antenna do not only affect the qubit 
itself, but also its readout circuitry. This means that 
the readout signal is usually covered by a wide spec- 
trum of resonances due to parasitic inductances and ca- 
pacitances of the DC SQUID leads. In order to avoid 
parasitic resonances, the DC SQUID bias resistors used 
for the resistive-bias measurements are fabricated par- 
tially on the sample chip (bias line: 250 Q. on-chip and 
1 kO off-chip; readout line: 2.25 kfi on-chip and Skfi off- 
chip). In this way, qubit, DC SQUID, and shunting ca- 
pacitor are located within a total length scale of about 
100 /im. Thus, most parasitic resonances in the relevant 
frequency range of a few gigahertz are strongly damped 
by the bias and readout resistors. The remaining reso- 
nances involve the shunting capacitor—, the inductance 
of the aluminum leads close to the DC SQUID, possi- 
ble box resonances, microscopic impurities in the sub- 
strate or the junctions, and, of course, the qubit. In the 
capacitive-bias detection scheme the small bias capacitor 
pushes parasitic resonances to higher frequencies. In ad- 
dition, smaller damping resistors are placed in a way that 
their low-frequency Johnson-Nyquist (thermal) noise is 
filtered out by the capacitors (bias line: 51117; readout 
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FIG. 4: (Color online) Filter characteristics of the DC SQUID 
bias line. The green (light gray) and the blue (dark gray) 
curves are recorded in the resistive and the capacitive bias 
configuration, respectively, (a) Shape of the bias current 
pulses. The switching pulse width above the hold level in the 
capacitive bias setup (~ 15 ns) is reduced by a factor of four as 
compared to the resistive bias setup (~ 60 ns) . (b) Transmis- 
sion spectra. In both cases, above 170 MHz the transmitted 
signal is below the noise floor of the network analyzer. In 
comparison to the resistive-bias setup, stronger attenuation 
is obtained below 15 MHz for the capacitive-bias setup. How- 
ever, the attenuation is weaker in the 15-170 MHz window. 
The width of this frequency range as well as the attenuation 
factor can be further optimized using a smaller bias capacitor. 



parasitic resonances. Secondly, although we are hmitcd 
to a 0.5 pF surface mount bias capacitor and although our 
qubit is energy relaxation limited at the optimal point, 
the relaxation time is not limited by the high-frequency 
bias current fluctuations (cf. Sec. IIVB|l . 

The capacitive setup permits us to obtain a narrow- 
bandwidth band-pass filter of 150 MHz width with a 
reasonably high center frequency of 90 MHz for the 
DC SQUID bias fine [cf. Fig.HJb)]. The minimum atten- 
uation is 55 dB. This means that we actually only open 
the frequency band which is absolutely necessary for the 
readout pulse. In future experiments, the use of smaller 
on-chip capacitors should allow us to design even nar- 
rower band-pass filters with lower minimum attenuation. 
We finally note that another advantage of the capaci- 
tive filter resides in the fact that it is not resistive and, 
as a consequence, especially suitable for mixing-chamber 
temperature applications. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

In this section, we present the experimental results of 
spectroscopy (Sec. IIV.^ . energy relaxation (Sec. IIVB|) . 
and dcphasing (Sec. IIV C[) measurements performed on 
our sample. In the discussion we first analyze in detail the 
results obtained with the capacitive-bias readout. Then, 
we briefly compare them to those achieved with the con- 
ventional resistive-bias method on the same qubit. In 
the case of the time domain measurements we focus our 
attention on the region near the degeneracy point of the 
qubit, where the coherence times are longest, and on the 
readout point. 

However, in the actual experiment the measurements 
with the well-known resistive-bias scheme are done flrst. 
From a technical point of view this has the advantage 
that the current-voltage characteristics of the readout 
DC SQUID can be obtained easily. The reason is that, 
as described in Sec. IIIIl only in the resistive-bias config- 
uration very low-frequency signals can be applied to the 
bias line. 



A. Switching discrimination and spectroscopy 



line: 1.5 kil). Because of technical reasons (realization of 
resistive and capacitive-bias configuration for one and the 
same qubit), in the capacitive-bias measurements only 
the shunting capacitance Csh is fabricated on-chip. For 
the bias and readout resistors / capacitors we use surface 
mount devices placed on the PCB as close as possible to 
the sample chip. In this case the on-chip resistors uti- 
lized in the resistive-bias measurements are shorted with 
approximately 1.5 mm long gold bonding wires. Fortu- 
nately, as it turns out, this docs not present a problem 
for our measurements. Firstly, the spatial enclosing of the 
sample to a length scale < 3 mm is still sufficient to avoid 



Before starting any experiments on the qubit dynam- 
ics, we need to confirm that the considerations made in 
Sec. mil regarding the feasibility of the capacitive bias- 
detection scheme are actually valid for the investigated 
three-junction flux qubit. In other words, we need to 
verify that the different voltage response signals of the 
DC SQUID, which are due to a switching or a non- 
switching event, can be clearly distinguished from each 
other. In our sample, this so-called switching discrim- 
ination is unambiguous provided that the length of the 
hold section of the bias pulse is at least 1.5 ^s (cf. Fig.E]). 
Additionally, we find that we do not detect any spurious 
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time 

FIG. 5: (Color online) Oscilloscope image of the readout 
DC SQUID signal (solid lines) vs. time for the integrated- 
pulse readout method. The switching rate is approximately 
50%. The broken red line marks the detection threshold 
used to discriminate between zero- voltage (no switching) and 
finite- voltage (switching) events. Note that no retrapping (re- 
turn to the zero-voltage state) occurs before the end of the 
1.5 fis hold pulse. The existence of two finite-voltage branches 
is due to a resonance in the current-voltage characteristic of 
the DC SQUID. The intermediate branch corresponds to the 
resonance step and the uppermost branch to the true super- 
conducting gap. 



2/c = 1202 ± 138 nA and, as a consequence^^, 

a = — ^ ^0.62 ±0.04. (3) 

2y^i-(/p//c)2 

This value is consistent with the result arcs — 0.63 ob- 
tained with the resistive-bias method and deviates only 
shghtly from the design value adesign = 0.7 (cf. Sec. 
Knowing a, we can calculate the Josephson coupling en- 
ergy E] /h = 299 ± 34 GHz. Here, the error is domi- 
nated by the uncertainty of the gap voltage due to sig- 
nificant fluctuations of the quasiparticle branch of the 
current- voltage charactristic of the DC SQUID detector. 
From previous measurements of junctions made with the 
same fabrication process^, the charging energy is esti- 
mated to be Ec/h ~ 6.4 ± 1.6 GHz, which gives a ratio 
E]/Ec = 47± 13. 

B. Energy relaxation 

In this section, we evaluate the results of the time 
domain measurements. We conduct energy relaxation, 
Rabi, Ramsey, and spin echo experiments using the pulse 



switching events on the oscilloscope when the readout 
pulse does not contain the switching segment. This con- 
firms that the detection of the qubit state occurs only 
during the short switching pulse. 

After this initial calibration we measure the qubit en- 
ergy diagram by means of pulsed microwave spectroscopy 
(cf. Fig. [H]). Each data point in the main body of Fig. [B] 
corresponds to the center frequency of the observed reso- 
nance peak. Let us first look at two typical spectroscopy 
traces which are displayed in the inset of Fig. [6l There, 
the qubit transition frequency is approximately 17 GHz. 
For the capacitive-bias method the qubit visibility, i.e., 
two times the height of the spectroscopy peak, is en- 
hanced as compared to the resistive-bias scheme. We find 
this increased visibility consistently in all spectroscopy 
and Rabi oscillation data, which has been recorded over 
a wide range of frequencies (data not shown) . More quan- 
titatively, the qubit visibility is 20-25% for the resistive- 
bias method and 28-36% for the capacitive-bias scheme. 
From the main body of Fig. [S] we find a tunnel coupling 
A/h = 4.22±0.01 GHz at the degeneracy point. Together 
with the latter result we can fit the fiux-dependence of the 
qubit resonance signal using a two-level system model. 
We include one-, two-, three-, and four-photon resonance 
peak positions into this fit and find a maximum circulat- 
ing current /p = 360 ± 1 nA. In order to estimate a, we 
make use of the fact that the two junctions of the read- 
out DC SQUID have the same layout as the two larger 
qubit junctions. In this way, we can estimate Ic from 
the measured normal resistance Rn = 258 f2 and gap 
voltage Vg = 390 ± 45 of the DC SQUID. Using the 
Ambegaokar-Baratoff relation 2/ci?n = we obtain 
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FIG. 6; (Color online) Qubit resonance frequency (various 
symbols) obtained from capacitive-bias spectroscopy at T ~ 
50 mK plotted vs. the external DC flux bias at T = 50 mK. 
Multi-photon resonances are observed up to the four-photon 
process. The solid lines are two-level system fits to the mea- 
sured data. Inset: Spectroscopy traces (switching probabil- 
ity vs. excitation frequency) for a qubit transition frequency 
i/Qi — 17 GHz. The visibility is two times the maximum am- 
plitude of the resonance peak. This is a typical example of the 
increased visibility of the capacitive-bias spectroscopy [blue 
(dark gray) curve] compared to the resistive-bias data [green 
(light gray) curve]. 
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FIG. 7: (Color online) Characteristic decay times of the qubit 
(left scale) plotted versus the external flux bias in the vicinity 
of the qubit degeneracy point. The data is obtained from time 
domain measurements with the capacitive-bias method. The 
magenta (middle gray) crosses represent the energy relaxation 
time Ti, the red (dark gray) open triangles the Rabi decay 
time TRabi, the green (light gray) open circles the Ramsey 
decay time T2R, and the dark blue (dark gray) open squares 
the spin echo decay time T2e- Also shown are the results 
of low-power spectroscopy (right scale). The light blue (light 
gray) plus signs represent the positions of the resonance peaks. 
The solid lines are fits to the data as discussed in the main 
text. In addition to the data shown in this plot, the results 
at the qubit readout point (5<l?x = 6.007 x 10"^ $0, foi = 
Eoi/h = 14.125 CHz) are used for the fits. 



sequences described in Sec. |TT] and Fig. [51 In our mea- 
surements we focus on the flux bias region around the 
degeneracy point, (5$x = ±6 x 10~'*$o corresponding 
to Eoi/h ~ 4 GHz, and on the readout point, S^^ = 
-6.007 X 10^3^0 corresponding to Eoi/h = 14.125 GHz. 
All results are summarized in Fig. [T] 

We now discuss in detail the energy relaxation of 
the flux qubit. According to Bloch-Redfleld theor y^^'^^ , 
which is applicable if the noise is short correlated and 
weak, the energy relaxation rate is given bj«2i 
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n /9e($x 



■ sm ( 



(4) 



Here, S^^{ujqi) = i ^S^^{~ujoi) + S^Hluoi)\ is the sym- 
metrized noise spectral density at the qubit transition 
frequency. The factor sin^ 6 = A^/[A^ + e(<i>x)^] can be 
understood intuitively because the Bloch-Rcdfield theory 
is based on a golden rule argument: Ti is certainly related 
to level transitions of the qubit and (1|(T^|0) = (Ola^ll) = 
sin 9 is the transition matrix element due to the interac- 
tion between qubit and noise. Also, the importance of 
the spectral density at the transition frequency becomes 
clear because only at this frequency state transitions can 



be induced. In order to calculate 5*1^(^01 ) from the flux 
noise spectral density S'<i,(a;oi), the latter has to be multi- 
plied by the square of the flux-to- frequency transfer func- 
tion C± = [L'(aa;/a$x)]_L = [h-^de{^^)/d<i>^]sm0 = 
{2Ip/h) sinO because the relaxation rate is determined 
only by the transverse fluctuations. Close to the degen- 
eracy point we have sinO ~ 1 and the relaxation rate is 
expected to show a very weak dependence on the exter- 
nally applied DC flux bias J^x- Looking at the experi- 
mental data of Fig. [71 we note that although there is some 
structure in the flux dependence of Ti, the relaxation 
times clearly exhibit a good agreement with the simple 
Bloch-Redfleld model near the degeneracy point. From 
a numerical flt we derive the flux noise spectral density 
to be Sl^iujoi ~ wa) = [(1.4±0.1) x 10-i°$o]^ Hz^\ 

With the help of Eq. ([4]) we can also estimate the re- 
laxation rate, which would be expected if the impedance 
Z{lu) of the DC SQUID biasing circuitry was the domi- 
nating source of energy relaxation. In order to do so, we 
have to use the voltage noise spectral density S'y (cjqi) as- 
sociated with the real part of the impedance, 3? {Z(cjoi)}. 
According to the fluctuation-dissipation theorem, for 
hujoi 3> k^T, which is true for our experimental situ- 
ation, we obtain iS'y(a;oi) = 2hujQi^{Z{LL!oi)} /2Tr. The 
flux noise density in the qubit loop due to Svi^Joi) is 
given by S'i>(woi) = M|q qi^i^T^S'y (wqi). Here, K'^ is 
the factor transforming the voltage noise spectral density 
into a current noise spectral density 5'/(a;oi) of the cur- 
rent circulating in the readout DC SQUID. The resulting 
S<t.{uJoi) for the qubit is simply obtained by multiplying 
Si{iOoi) by the quantity M|q which is the square of 
the mutual inductance between the readout DC SQUID 
and the qubit. Hence, we obtain 
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In the equation above, /gq = <I'sq/^'o = 0.97 x 1.5 is 
the magnetic frustration of the DC SQUID due to the 
flux threading its loop and Isq is the transport cur- 
rent through the DC SQUID. From the experimental 
parameters of the bias circuit (cf. Fig. [5]) we flrst cal- 
culate 3?{Z(woi)} (cf. Fig. [SD. Using Eqs. © and (O 
we can then estimate the lower limit for the relaxation 
time caused by the readout process. In this case, we 
can approximate /gq with the DC SQUID switching cur- 
rent /sw — 120 nA, obtaining Ti^z{ujoi=uia) ?i 2(US at the 
degeneracy point. This estimate is already much larger 
than the actually measured decay time Ti = 82 ns. Fur- 
thermore, in our Ti measurement protocol the bias cur- 
rent is non-zero only at the readout time (cf. Fig.[5|). For 
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FIG. 8: Equivalent circuit of the DC SQUID bias line. The ca- 
bles connecting the bias resistor to the voltage source are not 
considered here because the corresponding modes are strongly 
damped, (a) Applying Norton's theorem, the voltage-driven 
circuit displayed in Figs. [T] and |3] is transformed into an equiv- 
alent current-driven circuit. Istc is the effective source current. 
The DC SQUID is modeled by its Josephson inductance^ 
Lj — 2.6 nH. Ziine is calculated using Kirchhoff 's laws, (b) 
Details of the bias line impedance Zune- The inductances 
Li = 100 pH, La = 10 pH, L3 = 60 pH, and L4 = 45 pH 
are determined from FastHenry^ simulations. The bias 
impedance Zbias depends on the actual measurement setup, 
(c) Details of Zbias for resistive- (top) and capacitive-bias con- 
figuration (bottom). LRbias = 2.7nH and Lcbias = 1.5 nH are 
the self-inductances of the on-chip resistor and of the bonding 
wire used to short the on-chip resistor, respectively. 



this reason the readout process itself can only contribute 
to the Umited qubit visibility. The actual DC SQUID 
transport current during the Ti measurement sequence 
is given by asymmetry currents or other noise currents. 
This means that our experimental situation is charac- 
terized by the condition Isq <C /sw, which implies an 
even higher expected Ti time. Thus, we conclude that 
the noise generated by the bias circuitry is not the dom- 
inating high-frequency noise source in our experiments. 
We find a similar resulii^i for the resistive-bias configura- 
tion, where the measured relaxation time Ti = 140 ns 
is also small compared to the calculated lower bound 
Ti,z{uioi=i^&) ^ 5^s. The slight differences in the Ti 
values measured with the capacitive- and the resistive- 
bias readout arc mainly due to changes in the microwave 
setup, which occur during the switching between the dif- 
ferent experimental configurations. This process involves 
the warm-up of the cryostat, sample remounting, open- 
ing and closing of microwave connectors, and a new cool- 
down of the refrigerator. 

In conclusion, we find that the current noise from the 
bias circuitry is not the dominating relaxation source, 
neither for the resistive- nor for the capacitive-bias detec- 
tion scheme. Even changes in the high-frequency setup. 
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FIG. 9: (Color online) Real part of the impedance 
K {Z{Lj/2n)} of the DC SQUID bias line plotted vs. frequency 
for various configurations. The blue (dark gray) and the green 
(light gray) lines show the results for capacitive- and resistive- 
bias configurations, respectively. The actual configurations 
used in the experiments ("actual sample", thick lines) are 
shown in Fig. [8] The thin lines show suggestions for im- 
proved configurations the reduce K{Z(a;/27r)} in the relevant 
frequency range form 2 GHz to 10 GHz. 



which arc due to the switching between the two meth- 
ods, affect the Ti times by less than a factor of two. 
However, in the capacitive setup the measured relaxation 
rates are more than a factor of 25 larger than expected 
from Eq. ([5]). The actually dominating noise sources for 
our qubit are discussed in detail in Sec. IIVDI In future 
works, if these presently dominating noise sources can be 
reduced to an extent that the bias current noise becomes 
relevant, the capacitive-bias scheme allows to engineer 
the circuit impedance and the resulting spectral noise 
density in an advantageous way. As shown in Fig. [51 in 
the relevant frequency band from 2 GHz to 10 GHz, an 
impedance (and thus a spectral density), which is almost 
flat and approximately two to three orders of magnitude 
smaller than that of the present measurement setup, can 
be realized by a tenfold reduction of the bias capacitance. 
In contrast to that, for the resistive-bias readout in the 
same frequency band such a reduction cannot even be 
achieved by a twentyfold increase of i?Rbias- In other 
words, the bias line noise can be reduced significantly 
by applying the bias current via a capacitor instead of 
a resistor. Note that there are no major technical prob- 
lems hindering a tenfold reduction of the bias capacitor, 
whereas increasing -Rnbias by a factor of 20 can gener- 
ate significant problems due to heating effects in the bias 
circuit. 



C. Dephasing 



In order to determine the dephasing rate ry(<I>x) = 
we study the decay of the (Tj;.j,-evolution of the qubit 
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state vector on the Bloch sphere. This requires the use 
of muhi-pulse sequences such as Ramsey and spin echo 
sequence (cf. Sec. ITll and Fig. [2]). The reason for this is 
that the switching probabihty of the readout DC SQUID 
only reflects the (j^-component of the Bloch vector. The 
length of the 7r/2- and 7r-pulses, i.e., pulses which ro- 
tate the Bloch vector about an axis in the o'a;^j,-plane by 
the angles 7r/2 and tt, respectively, depends on the mi- 
crowave pulse amplitude and is determined experimen- 
tally from driven Rabi oscillations as described in Sec. |TT] 
(data not shown) . The azimuthal angle of the rotational 
axis in the (Ta;^y-plane is controlled by the relative phase 
of the pulses. In reality, the microwave control pulses 
are not perfect: The tipping angle can deviate from the 
desired value and the frequency can be detuned from the 
qubit resonance frequency. Neglecting relaxation and de- 
phasing, the Ramsey signal i?^^ and the spin echo signal 
^01,02 can be described by a series of rotations on the 
Bloch sphere: 

xR^^^s[Q.^)R^{2'K5t)Ro^s{n^/2) xq]^ 
and (7) 
= {5) 

= [R^,Xn^/2)Rz{2T:5t)RaXn^/2) xq]^ . 

Here, Xq = (0,0,1) is the Bloch vector describing the 
initial state of the qubit. The rotation matrices R^p^s 
and R, are those defined in Eq. They describe 

the Rabi dynamics and free evolution of the qubit state 

vector, respectively. ri7r/2 = f \/ ^ + {S/fR)'^ + ct* and 

il^r = t: \J~l~\~{sJu^ + P* are the angles by which the 
tt/2- and 7r-pulscs determined at resonance actually ro- 
tate the qubit Bloch vector in the presence of a detun- 
ing S and arbitrary pulse length imperfections a* and 
f3* . The phase of the initial 7r/2-pulse is chosen to be 
0°, that of the final 7r/2-pulsc 0i, and that of the in- 
termediate TT-pulse 02- In our experiments, we choose 
the microwave power such that the corresponding Rabi 
frequency is i>^^ ~ 65 MHz. This results in a 7r-pulse 
length = 2T7r/2 — 8 ns. A more detailed analysis of 
the measured Rabi oscillations shows that a* , P* < 5°. 
For perfect pulses {a* = (3* ~ 0°) with vanishing detun- 
ing (S = 0) and zero relative pulse phase {(pi = (j)2 = 0°), 
the measured signals should be Rq — —1 and Eq,q = +1 
independent of t. This means that in presence of decoher- 
ence we should directly observe the decay envelope, from 
which the decay time can be readily extracted. How- 
ever, in practice 6, a* , or /?* are non-zero, giving rise to 
time-independent (offsets) and time-dependent (oscilla- 
tions) changes in the signal. In some cases this makes 
the proper determination of the decay times rather diffi- 
cult. 



By adding or subtracting traces taken with differ- 
ent relative phases 0i and (j)2 the "ideal pulse" signals 
can be restored. This procedure is called phase-cycling 
method^!. In the case of the Ramsey traces, which 
are deliberately measured with detuning ~ 50 MHz, 
we can still cancel offsets by considering Rq — i?i8o or 
Rgo^ R270- There arc two important offsets: a small one 
(< 1%) due to pulse imperfections and the calibration 
offset of the switching probability. The latter is usu- 
ally chosen to be about 50%. In this way, the switching 
probability is also close to 50%, when the qubit is in a 
50/50- mixture of ground and excited state. In order to 
improve the signal-to-noise ratio we use all four available 
Ramsey traces 

R = R{S) = Ro — Riso + R90 — -R270 (8) 

for the data analysis (cf. Fig. [TO)). 

In the spin echo measurements there is no deliberate 
detuning {6 < 20 MHz) and we are able to keep only 
the non-oscillating spin echo term using the eight phase 
combinations 

E = E{S) 

= Eofl — Eiso^o ^ £'0,90 + £'l80,90 

+ ^'0480 — -£'180,180 — £'0,270 + ^'180, 270 7 (9) 

which of course decays with some (usually exponential) 
envelope in the presence of decohcrence. With the aid of 
numerical simulations we confirmed that, for the usual 
conditions of our spin echo experiments, these corrections 
are small. In fact, they disappear in the noise floor of 
the switching probability measurements. However, when 
looking at the corrected Ramsey and spin echo traces 
very close to the degeneracy point (cf. Fig. [Tl|), it be- 
comes evident that there are beatings which cannot be 
corrected by the phase-cycling method. These beatings 
arc discussed in more detail in Sec IVl 

We now turn to the analysis of the corrected Ram- 
sey and spin echo data, which provides information on 
the ^x,y- or transverse dynamics of the qubit and, as a 
consequence, on the spectral density of the phase noise 
affecting the qubit. The transverse signals have two main 
decay components. Firstly, there is the energy relaxation 
already discussed in detail in Sec. IIVBI which is obvi- 
ously still present. Secondly, there are low-frequency ran- 
dom phase fluctuations, which destroy the phase coher- 
ence of the qubit. When the energy relaxation is caused 
by regular high-frequency noise2ii2£, the total rates of the 
Ramsey (r2R) and spin echo (r2E) decay are given by 

Fi Fi 
r2R = — + F^R and F2E = + F^e , (10) 

where the contributions F^r and F^e are called the pure 
dephasing rates. The factor 1/2 is due to the fact that 
Fi is defined as an energy decay, whereas F2 and F;^ are 
defined as amplitude decays. Under the assumption of 
a regular ("white") noise spectral density Suj{w) at low 
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FIG. 10: (Color online) Ramsey decay traces measured at two 
different flux values very close to the degeneracy point using 
the phase-cycling method, r is the free evolution time. For 
the data [blue (dark grey) crosses] and the traces calculated 
according to Eqs. ((Tjl and ([8]) the same labeling is utilized. 
The solid green (light grey) lines are fits to the data using a 
split-peak model in combination with an exponentially decay- 
ing envelope. The observed decay times are (a) T2K = 75±4 ns 
and (b) T2R = 84 ± 5 ns. 



frequencies a; ~ 0, we can apply Bloch-Redfield theory 
agai n^^'^^ . In this way, we obtain a simple exponential 
decay envelope for the pure dcphasing and the decay rate 
is 



D 



s. 



cos 9 



S^^'icu = 0) . 



(11) 



In contrast to the energy relaxation, the flux-to- 
frequcncy transfer function C|| = [D{du}/d^x)] = 

[h-'^de{^y,)/d'P^]cos9 = {2Ip/h)cos9 has to be used 
since the dcphasing rate is determined by the longitudi- 
nal fluctuations. The term — (Ojo-xlO) = (l|o-z|l) = cos6' 
can also be understood intuitively because no level tran- 



FIG. 11: (Color online) Spin echo decay traces measured 
very close to the degeneracy point using the phase-cycling 
method. The same labeling is used for the data and for 
the traces calculated according to Eqs. Q and (a) 
The aperiodic low-frequency beatings (< 2 MHz) are inde- 
pendent of the relative phase of the 7r-pulse and do not can- 
cel in the corrected trace [thick red (dark gray) line]. Sim- 
ilar to the Ramsey trace R, the spin echo difference traces 
-E^j = (i?0i,i8o — E^-^fi) (thin lines) have a built-in offset 
correction, (b) Corrected spin echo decay obtained from (a). 
The solid magenta (medium gray) line is a fit to a simple ex- 
ponential envelope function. From all data points [blue (dark 
gray) crosses] only those with absolute values close to the en- 
velope [green (light gray) plus signs] are used for the fit. The 
amplitude is treated as a fixed parameter which is evaluated 
from the Ramsey amplitude and the Rabi decay time (Rabi 
decay data not shown). 



sitions are induced by the low-frequency phase noise. 
Furthermore, we notice that the flux dependence of the 
dcphasing rate of the qubit is dominated by the factor 
cos^ 6* = e((f>x)^/[A^ -I- e(<i>x)^] , which is proportional 
to ((5<i>x)^ close to the degeneracy point and approaches 
unity far away from it. We also notice that the dcphas- 
ing rate is determined by S^^{uj = 0). This means that 
l//-noise is of particular importance for the phase co- 
herence of the qubit. However, l//-noise obviously vi- 
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olates the assumption of regularity and one has to cal- 
culate the actual noise spectral density from the total 
accumulated phase^i. The result depends on the ap- 
plied control sequence (Ramsey or spin echo), the order 
of the coupling (linear or quadratic), and the choice of 
the frequency cutoff (sharp or crossover to 1//^). It can 
be Gaussian, exponential of (where a = 3,4), or 
even algebraic. Unfortunately, the limited overall data 
quality of the traces shown in Fig. [TT] (short coherence 
times, switching probability fluctuations, and beatings) 
does not allow a detailed trace-by-trace analysis of the 
experimental data. Consequently, we follow a different 
strategy. We determine the decay rate r2 from a simple 
exponential fit of the trace envelope. Because of the ape- 
riodic low-frequency beatings we fit the spin echo traces 
taking the amplitude as a fixed parameter, which is deter- 
mined from the Ramsey amplitude and the Rabi decay 
time. Then, wc fit the flux dependence of r2 with the 
expression 



r2(<5$x) 



ri((5$x) 



+ r^^(5<i'x) 



(12) 



Here, the pure dephasing rate has a Bloch-Redficld 
(white flux noise) contribution F^^, a Gaussian (1// flux 
1 / f 

noise) contribution F^' , and a constant (quadratic cou- 
pling, flux-independent sources) contribution F|^. For a 
1/ /-noise spectral density 



Sl/^{u)^A/\Lo\ 



one obtains21 



9$. 



cos 9 



V^ln2 



(13) 



(14) 



Strictly speaking, Eq. p4|) is only true for the spin echo 
decay. However, the correction necessary for the Ramsey 
decay rate is only of logarithmic order. Its effect is to 
increase the impact of the noise by a logarithmic factor, 
which can be regarded as a constant in our considera- 
tions. This means that the T'2R-fit produces a Ramsey 
l//-noise magnitude which is by itself unphysical. 
Nevertheless, is useful to understand the filtering ef- 
fects of the spin echo sequence and estimate the infrared 
cutoff frequency wir of the l//-noise spectrum. 

In our experiments, the Ramsey and spin echo data is 
taken in the vicinity of the degeneracy point ((5<I>x = ±4x 
10"''<I>o, i^ai = 4.22 - 4.30 GHz) and at the qubit read- 
out point ((5$x = -6.007 x IQ-^^g, i,^^ = 14.125 GHz). 
The details of the fitting procedure arc different for Ram- 
sey and spin echo data. The Ramsey traces can be fit- 
ted properly assuming a simple model in which the main 
qubit resonance is split into two peaks. Both peaks have 
equal amplitudes and are detuned by the quantities Si 
and 62 from the main resonance. Using R{6) of Eq. ([8]), 
the total decay function consists of beating oscillations 
Rat = [^(^i) + -R(<52)]/2 multiplied by a simple expo- 
nential decay envelope. The observed peak splittings 



1^1 — (52 1 vary between 3 and 15 MHz around the degener- 
acy point. At the readout point no splitting is observed. 
In contrast to the Ramsey beatings, the spin echo beating 
cannot be fitted using a phenomenological model. From 
Fig. [TlT a) one can see that only the smaller wiggles are 
canceled by the phase-cycling method. The long-time- 
scale beatings coincide for all pulse sequences regardless 
of the relative phase shifts. For this reason these beat- 
ings are not considered as an artifact due to imperfect 
pulses. Consequently, we fit a simple exponential decay 
to the envelope of the long-time tail [green (light gray) 
plus symbols in Fig. [TlTb)] of the total spin echo trace 
[blue (dark gray) crosses in Fig. [TlT b)]. A substantial 
part of the error bars of T2E comes from the fact that it 
is not always easy to define this tail with high accuracy. 

The flux dependence of the decay times T2E = F^g 
and r2R = r2R displayed in Fig. [71 Obviously, 
the spin echo decay time is considerably longer than 
the Ramsey decay time. This indicates the presence 
of low-frequency noise in the system, which is partially 
canceled by the echo sequence. Fitting the T2E data 
with Eq. p2)) . we obtain a l//-flux-noise magnitude 

A = [(4.3 ± 0.7) X 10"*^$o]^- Right at the degeneracy 

point we find Tse 2Ti and = (FOe)"' ~ 2 fis, 
i.e., the total qubit coherence is limited by energy re- 
laxation. These results coincide with those obtained 
with the resistive-bias readout method^!, although the 
maximum r2E is by a factor of two larger in the latter 
case. This is not surprising, considering the Ti limitation 
and the fact that already the Ti analysis (cf. Sec. IIVB[) 
suggests an increase of high-frequency noise due to the 
switching from the resistive- to the capacitive-readout 
scheme. Furthermore, we obtain a finite Bloch-Redficld 
contribution S'|1(w = 0) = [(2.1±0.1) x 10-i°$o] ^Hz"^ 
to the spectral flux noise density. In the case of the Ram- 
sey decay, we also find that the flux dependence of the de- 
cay time T2R is consistent with the existence of l//-noise 

(Ar = [(2.0 ± 0.2) X lO-^'fo]') and T^^ {Vl^y' 
200 ns. However, wc did not find a significant Bloch- 
Redfield contribution 5'f'j^(cj = 0) in the entire fiux in- 
terval from (5$x = -6.007 x 10"^$o to +4 x 10"*$o- In 
contrast, for the spin echo decay a significant part of the 
low- frequency l//-noise is canceled (Ar/A ~ 5) by the 
refocusing effect of the intermediate 7r-pulse. 



D. Energy relaxation vs. dephasing and noise 
sources 



After the quantitative analysis of the influence of the 
external magnetic flux noise on the qubit decoherence 
performed in Sees. IIV Bl and llVCi we are now ready to 
put together the results into a general picture in this sec- 
tion. Let us recall that there are three classes of decay 
times analyzed in this work. The first one consists of the 
so-called Bloch-Redfield terms (Fi, F^^), which can be 
derived with a golden rule type of argument. The sec- 
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FIG. 12: (Color online) Ramsey and spin echo filtering func- 
tions. 

1 / f 

ond is due to the impact of 1 //-noise {TiJ ) and, finally, 
there is a contribution F^, which does not depend on 
the external flux. In order to understand the interplay 
of these three terms we have to consider their derivation 
in more detail. Focusing on the dephasing first, under 
the assumption of a Gaussian distribution of the noise 
amplitudes we obtain the decay envelope^ii^iiH 



fir) 



9$^ 



cos 6 



X / S^{u;)F{u;T/2)dLj } . (15) 



The filtering functions F ; 
spin echo sequences, 



FR,Ei'^) for the Ramsey and 



sin^w 
t'B.[uj) = and 



Fe(w) 



sin"* {Q/2) 



(16) 
(17) 



arc plotted in Fig. [T^l Here, we use uj = ujt/2. We 
can see that the Ramsey sequence is most sensitive to 
noise close to a; = 0. On the contrary, the spin echo se- 
quence is completely insensitive to zero-frequency noise. 
However, the corresponding filtering function exhibits a 
maximum at the finite frequency u ~ 4/r. Noise at this 
frequency just cancels the refocusing effect of the inter- 
mediate TT-pulse. In other words, l//-noise is strongly 
reduced by the spin echo sequence, while noise which is 
white in the angular frequency band between zero and 
Uc ~ 12/r remains unaffected. For typical time scales of 
our experiments, r « Ti,T2 ~ 100ns we find that the 
critical angular frequency ujc is still much smaller than 
the angular frequency uja corresponding to the qubit gap 
A. 

Both the l//-noise decay rate of Eq. (fTil) and the 
Gaussian decay law can be calculated straightforwardly 



from Eqs. and The Bloch-Redficld de- 

phasing rate of Eq. (fTT|) has to be treated carefully: The 
golden rule argument is only valid for times long com- 
pared to uj^^. Since the dephasing is dominated by the 
low-frequency contribution (cu — > 0) we cannot apply this 
argument to our time traces (typically 1 — 450 ns long) . 
Instead, we notice that Eq. (|lip and a simple exponen- 
tial decay law can be derived rigorously from Eqs. p5p 
and (fTT)) assuming a frequency-independent spectral den- 
sity. Thus, a finite F^^ contribution indicates the exis- 
tence of white noise. If both white and 1//- noise are 
present, near the degeneracy point the peaked shape 



of F 



cx (|cos0|) dominates over the rounded 

shape of (F™) ^ cx (| 008^6*1) ^. The peaked shape of 
T2K and T2E in Fig. [7] thus represents a clear evidence of 
the influence of l//-noise on our qubit. Near the read- 
out point = 6.007 x 10^^<I>o, we have \cos9\ ~ cos^ 9 
and the influence of the white noise on T2E becomes com- 
parable to that of the 1//- noise. We find that it is not 
possible to properly fit the T2E-data (including the read- 
out point), when omitting either the Bloch-Redficld or 
the l//-term. This means that wc find a transition from 
a 1/ /-noise dominated regime in the direct vicinity of the 
degeneracy point to a another regime near the readout 
point, where the influence of 1//- and white noise is com- 
parable. In the Ramsey data the l//-noise contribution 
is clearly dominant because of the singularity at a; = 0. 
For this reason it is difficult to extract the small white 
noise contribution from the data with a meaningful error 
bar. The fact that the spin echo sequence strongly re- 
duces l//-noise, but does not affect white noise, allows 
us to detect the latter. 

In order to calculate the true Ramsey decay law due 
to 1 //-noise from Eq. ([T5)) . we need to assume at least 
an infrared cutoff tJiR. The corresponding decay time is 
then 



5e($x) 



cos 



Aln- 



1 



Wirt 



(18) 



Gonsidering a maximal time trace length of about tq 
300 ns, we obtain 



2ttto exp 



In 2 



0.2 Hz. 



(19) 



This result has the same order of magnitude as the in- 
verse time it takes to record a single averaged data point. 
However, the exponential in Eq. (fTOl) makes the result 
exponentially sensitive to measurement errors. Consid- 
ering that a better time scale for (a;iR/27r) ^ is the time 
to record an entire decay trace (~ 1 h) , we probably over- 
estimate the infrared cutoff frequency slightly. 

Furthermore, we want to point out that the white noise 
spectral density deduced from the dephasing, S'f^, is 
of the same order as the Bloch-Redficld spectral den- 
sity S'^^(wa) inferred from the energy relaxation rates. 
The energy relaxation rate of Eq. (jj]) is calculated with 
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the golden rule assuming weak noise characterized by 
a smooth spectral density at o^a- In contrast to the 
Bloch-Redfield dephasing rate, the long time condition 
t » (w^^ ^ 250 ps) is fulfilled for our experimental time 
scales of 1 — 300 ns. The fact that we find S'f^ ~ 

Sij>^{uj/s,) ^ (wa) suggests that the white noise ob- 
served in the T2E-analysis is also limiting the energy re- 
laxation times. Considering that the noise created by the 
bias circuit cannot explain the observed relaxation rates 
(cf . Sec. IIVB|) and that Ti is quite sensitive to changes in 
the microwave setup, the natural candidate for the dom- 
inating white noise source is the high-frequency line. In 
contrast to that, l//-noisc can be modeled with an en- 
semble of microscopic fluctuators close to or within the 
qubit junctions^. In Sec. |Vl we show that we can ob- 
serve the effect of at least one such fluctuator in our time 
domain traces. 

Let us now compare the performance of our qubit to 
that of flux qubits with a galvanic coupling to the readout 
DC SQUID, which have been measured recently^^i^. As 
a prototypical ex amp le, we compare our results to those 
presented in Ref. |26| . The most striking similarity is the 
relaxation limitation at the optimal point. Furthermore, 
the mutual inductance between the qubit and the high- 
frequency line as well as the bias line filtering are almost 
the same. However, the Ti-time of our qubit is about one 
order of magnitude smaller. These observations immedi- 
ately imply that our qubit is more susceptible to high- 
frequency noise. We find two main reasons to explain this 
fact. First, any noise of given amplitude couples more 
strongly to our qubit because the persistent current is 
about twice as large as the one in Ref. [2^. Second, and 
more important, the total coupling between qubit and 
DC SQUID is much lower in our design (no kinetic mu- 
tual inductance and smaller DC SQUID junctions). For 
manipulations in the vicinity of the optimal point we ap- 
ply an adiabatic shift pulse (cf. Scc lII C[) during the qubit 
operation via the virtually unfiltcrcd high-frequency line. 
In contrast, in Ref. [2^ the readout pulse, which was sent 
through the strongly low-pass filtered bias line, causes a 
flux shift due to the large coupling. This shift does not 
take place during the qubit operation. In this respect the 
galvanic coupling method turns out to be advantageous. 
On the other hand, the controlled application of the shift 
pulse offers a large flexibility in the experiments. In par- 
ticular, the possibility to operate the superconducting 
magnet in the persistent mode is attractive. For this rea- 
son, we plan to use an additional low-pass filtered control 
line for the shift pulse in future experiments. 

Finally, other possible noise sources that can deteri- 
orate the qubit relaxation time are dielectric losses in 
the oxide barrier of the junctions^. In a three-junction 
qubit the whole loop acts as a very big fourth junction. 
Thus, the controlled introduction of a small fourth junc- 
tion as well as the use of small loop sizes (as it was done 
in Refs. [2^ and [4ll ) might be advantageous because the 
effective amount of dielectric material in the loop is sig- 
nificantly reduced. 



V. SPIN ECHO AND RAMSEY BEATINGS 



In this section, we analyze in more detail the beat- 
ings which we observe in the Ramsey and spin echo data 
shown in Figs. [TO] and [TT] In the previous sections, we 
already excluded trivial origins of such beatings. Firstly, 
as pointed out in Sec. llVCi they cannot be caused by 
imperfect control pulses because we applied the phase- 
cycling method. Also, because of the spatial enclosing of 
the sample in combination with the filtering conditions 
of the measurement lines, parasitic resonances in the ex- 
perimental setup present only a remote possibility (cf. 
Sec. lIIIj) . Finally, as we can see in Fig.l^l the electromag- 
netic environment in the immediate vicinity of the qubit 
does not exhibit any prominent mode with a frequency 
close to the qubit gap frequency A/h. In the following, 
we suggest a different explanation for the origin of the 
observed beatings. This explanation exploits the formal- 
ism of cavity quantum electrodynamics, which, in gen- 
eral, describes the interaction between atoms and light. 
In quantum optics, the coupling between an atom and 
the photonic mode of a cavity can give rise to beating 
signals^i^. In our case, the qubit acts as an artificial 
atom and the cavity is represented by a quantum har- 
monic oscillator located on the sample chip close to the 
qubit. As we will see later, the nature of this resonator 
is microscopic in the sense that its coupling to the mi- 
crowave drive is much smaller than the qubit-microwave 
coupling. However, the coupling strength between qubit 
and resonator is larger than the characteristic decay rates 
each of them. Under this assumption the system can no 
longer be described in terms of a spin-boson model. In- 
stead, the qubit-oscillator system has to be treated as a 
single quantum system. All our arguments given above 
are supported by numerical simulations, which qualita- 
tively reproduce the main experimental features. More- 
over, it is notable that our discussion is not limited to 
a specific type of qubit or readout method. In addition, 
when the oscillator is undriven and its thermal popula- 
tion negligible, it qualitatively behaves very similarly to 
a two-level system. This statement agrees well with our 
simulations. For completeness, we state that the effect 
of two-level fluctuators was also addressed in a previous 
work^, but with a different focus. 



The assumption of a harmonic oscillator interacting 
with the qubit is motivated by spectroscopy experiments 
performed with the resistive-bias readout. There, we 
observe a series of peaks which are reminiscent of the 
two-photon blue sideband transition signal of a qubit- 
resonator system, in which the harmonic oscillator has a 
frequency of — 4 GHz ~ A//i. These peaks have the 
same flux-dependence as the normal two-photon peaks, 
however with a frequency offset (cf . inset of Fig. [Ti|) . To 
model this scenario properly, we make use of the follow- 
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FIG. 13: (Color online) Simulated periodic Ramsey beatings 
due to a qubit-resonator interaction. The resonator (a quan- 
tum harmonic oscillator) is truncated to 31 bases. The sys- 
tem parameters are chosen to be similar to those found in 
the experiment: The qubit persistent current 7p is 370 nA, 
the qubit gap frequency A/h is 4 GHz, and the resonator 
frequency Vr is 4 GHz (resonance condition). The maxi- 
mum free evolution time r is 450 ns. The microwave excita- 
tion frequency is slightly detuned from that of the oscillator 
{v = i^Ramsoy = 4.05 GHz). The coupling constants are chosen 
in a way that the results resemble our experimental data, (a) 
Variation of the Ramsey signal with increasing gmw.r. The 
other parameters are flqb.r = 2 MHz and (?mw,qb = 80 MHz 
for the displayed traces, (b) Variation of the Ramsey signal 
with (?mw,r: The beating frequency increases with Pqb,r. Our 
experimental results indicate gqb,r — 2-6 MHz. The other pa- 
rameters are gmw.r ~ and ^mw.qb ~ 80 MHz. 



ing Hamiltonian, which we write in the qubit eigenbasis: 



+ i<7qb,r (S,^ — a) (cos 6 (Jz ~ sin 6 Ox ) 
— 5mw,qb COS (^Txvt) (cos 9 Gz — sin 6 Ox ) 

+ iffmw.r COS {2'ITVt) (o) — o) . (20) 

Here, the first two lines denote an undriven qubit- 
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FIG. 14: (Color online) Simulated aperiodic spin echo beat- 
ings due to the interaction of the qubit with a single harmonic 
oscillator. The total free evolution time is r and pmw.r ~ 0. 
Unless otherwise stated, the system parameters are the same 
as in Fig. 1131 The microwave driving is resonant with the os- 
cillator {v = !^Echo = 4 GHz). Depending on (?qb,r the beating 
spin echo trace crosses or not the zero axis within the cho- 
sen time window. If, within the frequency range of interest, 
the qubit is coupled to several (but not too many) spurious 
resonators, their influence has a complex dependence on the 
external flux bias. For this reason, the presence or absence of 
the Ramsey and spin echo beatings does not posses any ap- 
parent flux bias dependence. However, for a certain flux bias 
the same beating pattern occurs reproducibly. Such a behav- 
ior agrees with our experimental observations (cf. Figs. llOl and 
lll|l . Inset: Resistive-bias spectroscopy close to the degener- 
acy point. The green (light gray) plus symbols correspond to 
the qubit resonance peaks. The blue (dark gray) crosses are 
consistent with the two-photon process of the blue sideband 
transition due to the interaction of the qubit with one specific 
parasitic resonator. 



resonator system, which is similar to the standard 
Jaynes-Cummings model. The last two lines describe a 
microwave driving at the frequency i^, which is different 
for the Ramsey {u = i^Ramscy) and spin echo {u = i^Echo) 
experiments. These driving terms are only present within 
the duration of the control pulses, and they are switched 
off during the free evolution periods. The operators 
and a are the usual boson creation and annihilation op- 
erators for a harmonic oscillator, and gqb.r, 5mw,qbj and 
5mw,r the coupling constants. The (cos CTz — sin ct^)- 
terms appear because of the rotation of the interaction 
Hamiltonian into the qubit eigenbasis. As a matter of 
fact, the (cos 9 az)-teini constitutes the main difference to 
the Jaynes-Cummings approximation and is accurately 
accounted for in our simulations. Since we are investi- 
gating a flux qubit. all interactions are proportional to 
pairwise products of the qubit current operator Iqh , the 
oscillator current operator /j- , and the classical microwave 
driving current cos (2-Ki't). Also, we can restrict our 
discussion to the case of pure linear coupling because the 
microwave power used in our experiments is relatively 
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small. 

Based on Eq. (PO)) . we simulate the Ramsey and spin 
echo sequences at the degeneracy point. We choose simi- 
lar parameters to those found in the experiments, but ne- 
glect the effects of decoherence for simplicity. The results 
are shown in Figs.[Tn]and[T31 Already at a first glance it is 
evident that the simulations can reproduce the two main 
features of the experimental data: There can be peri- 
odic Ramsey beatings [Fig. fTSl b)] and, at the same time, 
there can be quite aperiodic spin echo beatings (Fig[T3]). 
The former are due to an on-resonance qubit-resonator 
interaction, whereas the latter are caused by the the re- 
focusing TT-pulsc, which generates an unwanted evolution 
in the resonator part of iJqb.r- When the resonator driv- 
ing strength guiw,T is increased from zero, the Ramsey 
beatings are washed out [Fig. [T^a)]. This indicates that 
in our experiments the resonator is a point-like micro- 
scopic entity (e.g., a magnetic impurity or a dangling 
bond in one of the qubit junctions), rather than an ex- 
tended LC-circuili^. We use the term "point-like" to 
denote that the resonator must have a small (compared 
to the qubit) effective area interacting with the exter- 
nal microwave drive. This means that the interaction 
constant (?mw,r ^ ffmw,qb, ffqb.r- Consequently, we choose 
.9mw,r = in the subsequent simulations [cf. Figs. [TMb) 
and [14]. Nevertheless, the spurious resonator is found 
to be stable enough to survive at least one thermal cy- 
cle. In fact, Ramsey beatings can also be observed in the 
resistive readout measurements (time domain data not 
shown, for spectroscopy data cf. inset of Fig. [M]). On 
closer inspection, we find that the resonance condition 
is very important for our simulation results. This might 
seem controversial, given that we observe the beatings 
over a considerable frequency range around the degener- 
acy point. However, we have to recall that we neglect 
decoherence here. As a consequence, in the simulations 
we have much sharper resonance peaks than in the ex- 
periments. Furthermore, we only account for a single os- 
cillator. In reality, there are probably several oscillators 
with slightly different frequencies close to 4 GHz and with 
different coupling constants. It is noteworthy to mention 
that there is no fundamental restriction for these oscil- 
lators to be present only in the frequency range around 
4 GHz. However, only there the qubit coherence time is 
long enough to allow for a clear observation of beatings 
or splittings. The parasitic resonators are also a possi- 
ble explanation for the fact that the measured relaxation 
times are much shorter than theoretically expected from 
the bias current fluctuations (cf. Sec. IIVB]) . 

The fact that the simulations suggest a point-like struc- 
ture as source of the beatings makes it difficult to de- 
cide whether we have to assume a harmonic oscillator 
or a two-level system. Indeed, in the simulations an 
undriven resonator behaves very similar to a two-level 
system. Nevertheless, we continue to use the term "res- 
onator" in order to keep the language simple and easy to 
understand. We also want to point out that in the ex- 
periments performed with the resistive-bias readout we 



find similar Ramsey splittings, but no clear spin echo 
beatings (data not shown). We attribute this to the fact 
that in the resistive experiments the visibility is worse, 
while the coherence times are longer. For the rather ape- 
riodic low- frequency spin echo beatings (< 2 MHz), this 
results in the tendency to observe the decay envelope 
only, whereas for the faster, periodic Ramsey beatings 
(~ 10 MHz) the beatings are still visible. Also, the res- 
onance condition discussed above is important for the 
observability of beatings. The reason is that different 
oscillators interact with the qubit via different coupling 
constants. However, beatings cannot be observed unless 
the resonator frequency falls into the frequency window 
in which the qubit coherence times are longer than the 
inverse of the coupling constant. This window is set by 
the qubit gap, which, due to thermal cycling, shifts by 
almost 10% from Apbias/^ = 3.9GHa21 in the resistive- 
bias experiments to Acbias/^ = A//i = 4.22 GHz in the 
capacitive-bias experiments. 



VI. CONCLUSIONS 

In summary, we have measured the coherent dynam- 
ics of a three- Josephson-junction superconducting flux 
qubit, which is characterized by a gap frequency A/h = 
4.22 GHz. We used the capacitive-bias detection method, 
a variant of the switching DC SQUID readout. This 
method provides a built-in band-pass filtering of the elec- 
tromagnetic environment of the qubit, thereby allow- 
ing for faster readout pulses and a reduction of low- 
frequency bias current noise. We use the adiabatic shift- 
pulse method to probe the qubit close to the degeneracy 
point and the phase-cycling technique to cancel the ef- 
fects of control pulse imperfections. Right at the degen- 
eracy point the coherence time of our qubit is limited by 
energy relaxation and we find Ti = 82 ns and ~ 2 /zs. 
The flux noise spectral density at the qubit transition fre- 
quency is SI^{u;a) = [(1.4 ±0.1) x 10-io$o] ^ Hz'^ In 
the close vicinity of the degeneracy point, the measured 
spin echo times exhibit a flux dependence that is expected 
for the case where flux noise with a l//-frequency spec- 
trum and a mag nitude A = [(4.3 ±0.7) x lO-^^o] is the 
dominant source of dephasing. Near the readout point, 
which in our case means far away from the degeneracy 
point, we also find a significant white noise contribution 
S^^ = [(2.1 ±0.1) X 10-i"$o]^Hz~\ From the com- 
parison between Ramsey and spin echo decay times we 
estimate an infrared cutoff wiR/27r ~ 0.2 Hz of the 1//- 
noise spectrum. 

On one and the same qubit, we compare the results 
obtained with the capacitive-bias readout to data taken 
with the conventional resistive-bias method. We find that 
the capacitivc data shows an increased visibility. On the 
other hand, a slightly reduced (less than a factor of two) 
Ti-time is observed. The 1/ /-noise magnitude deduced 

from the spin echo decay is about (lO~®<f>o)^ in both 
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FIG. 15: (Color online) Sketch of the relevant noise sources af- 
fecting the qubit: We find white noise (arrows), 1//- noise [en- 
semble of blue (dark gray) dots] , and a few fluctuators coupled 
especially strongly to the qubit [blue (dark gray) parabolas]. 
In our experimental scenario, the high-frequency bias noise 
of the DC SQUID [green (light gray) arrow] is much smaller 
than the noise generated by the high-frequency line [magenta 
(medium gray) arrow]. 

setups (cf. Ref. [13 for more details on the resistive-bias 
results). This value was also found in a study where 
the flux qubit had a galvanic connection to the readout 
DC SQUID2^. Together with the detailed study on the 
influence of high-frequency environmental bias noise on 
the energy relaxation time Ti, our work clearly shows 
that external noise generated by the DC SQUID bias 
circuitry is not the dominant noise source for either T2- 
or Ti-time. We attribute the difference in the Ti-times 
obtained with the capacitive- and resistive-bias measure- 
ments to changes in the high-frequency setup, which have 
occurred during the experiments in the process of switch- 
ing between the two methods. Finally, we flnd that in the 
frequency range of 2 — 10 GHz, which is most important 
for flux qubit operation, the capacitive-bias scheme al- 
lows one to easily engineer an environment with a low 
and almost flat noise spectral density of the bias current 
fluctuations. This is not possible for the resistive-bias 
readout. In other words, if, in the future, the presently 
dominating noise sources can be reduced to a level that 
bias circuit noise becomes the limiting factor, the elec- 
tromagnetic environment of the capacitive-bias method 
is advantageous. 

Near the qubit degeneracy point, we find beatings both 
in the Ramsey and the spin echo traces. These beatings 



constitute an experimental evidence of the interaction of 
the qubit with one (or a few) harmonic oscillators or two- 
level fluctuators, which are approximately on resonance 
with the qubit. We support this argument with the aid 
of numerical simulations, which qualitatively explain our 
data. In particular, these simulations indicate that the 
perturbing oscillators are point-like, i.e., their effective 
area, and thus their coupling strength to the microwave 
driving, is small compared to the qubit loop area. 

Finally, comparing the white-noise contribution to the 
dephasing with the l//-contribution and with the noise 
responsible for the energy relaxation, we flnd ~ 

'S'^^(wa) ^ Slf'^iLu^). This means that, according to 
our experimental results, there are presently two main 
noise sources limiting the qubit coherence (cf. Fig llSp . 
The first one is white noise, which can exist up to the 
qubit gap angular frequency wa- This noise represents 
the main reason for the relaxation limitation at the de- 
generacy point and also gives a significant contribution 
to the dephasing of the spin echo signal far away from 
the degeneracy point. This noise is not generated by 
the DC SQUID bias circuit and, as a consequence, most 
probably couples via the high-frequency line to the qubit. 
In contrast to that, we observe the effect of l//-noise on 
the qubit dephasing in the close vicinity of the degen- 
eracy point. Because of the observed beatings in the 
Ramsey and spin echo signals, we attribute this noise to 
an ensemble of microscopic fiuctuators or resonators lo- 
cated on the sample chip close to or even within the qubit 
junctions. 
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